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\& • We discuss a new mechanism which can solve the crisis in standard big- 
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r~| ■ bang nucleosynthesis. A longdived particle X (10 4 sec ^$ tx ^ 10 6 sec) which 

a.: 

q^ 1 decays into photon(s) will induce cascade photons, and destroy significant 

<D : 

• amounts of D and 3 He without destroying 4 He or too much 7 Li. We numer- 
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ically investigate this process and derive a constraint on the properties of X 
such that the theoretical values of the primordial light-element abundances 
agree with observation. We also present some candidates for the unstable 
particle X. 
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The abundances of the light elements are a crucial test of the big-bang nucleosynthesis 
(BBN) scenario, since BBN precisely predicts the primordial (i.e., before stellar processing) 
abundances of D, 3 He, 4 He, and 7 Li ||T|. In the past, improvements in the theoretical and 
observational estimates of the light elements seemed to give us a better confirmation of the 
standard BBN scenario. 

Recently, however, observations seem to conflict with the theoretical predictions of stan- 
dard BBN: if we determine the baryon-to-photon ratio rj by using the observed 4 He abun- 
dance, then the primordial abundances of D and 3 He predicted from the standard BBN 
calculation are too large. In particular, based on a statistical argument, Hata et al. claimed 
that the standard BBN scenario is inconsistent with observations at a very high confidence 
level [||. Furthermore, Copi et al. pointed out that the systematic error in the 4 He observa- 
tion has likely been underestimated if the standard BBN scenario is correct. These arguments 
suggest the need for some modification of the theoretical or observational estimations of the 
primordial abundances. 

At present, there are several attempts to solve the discrepancy, though none of them 
has been confirmed. One way is to adopt some modification of the observed abundance. In 
particular, if we assume an extra systematic error in the estimation of the 4 He abundance 
or if we trust the higher measurements of D in the primordial Hi cloud, then we will have 
better agreement between theory and observation || . Alternatively, the models of chemical 
and stellar evolution of the light elements may need some change. The other way to solve 
the problem is to modify the standard scenario of BBN. For example, one can reduce the 4 He 
abundance by allowing degenerate electron neutrinos |4j], or by assuming a massive unstable 
neutrino with mass about 1 MeV and lifetime about 1 sec ||. Since BBN, together with the 
Hubble expansion of the universe and the cosmic microwave background radiation (CMBR), 
is one of the most essential predictions of the big-bang cosmology, it is important to clarify 
the origin of the discrepancy. 

In this letter, we propose another new mechanism to solve the crisis by reducing the 
predicted D and 3 He abundance. Essentially, if some exotic particle X decays into photons 

1 



well after BBN, then the predictions of standard BBN can be modified. Conventionally, such 
radiatively decaying particles have been regarded as sources of cosmological difficulties . 
However, we can actually have a better scenario of BBN if the cascade photons can destroy 
significant amounts of D and 3 He without affecting 4 He and 7 Li. Therefore, in this letter, 
we derive a constraint on the properties of X which makes the BBN scenario viable. We 
also suggest several candidates for the unstable particle X based on realistic particle physics 
models. 

First of all, we briefly review the effects of the radiative decay of X. Once a high 
energy photon is emitted into the thermal bath, cascade processes are induced. The high 
energy photon, electron, and positron spectra are formed accordingly. Cascade photons then 
induce photodissociation of the light elements. Thus, if the radiative decay of X induces 
many cascade photons at a time t after BBN starts (t ^ lsec), then the predictions of 
standard BBN can be modified. In the photon cascade process, double photon pair creation 
(7 + 7bg — * e + + e~, with 7 B q being the photon in thermal background) plays an important 
role: almost all the photons with energy larger than e t h r ~ m 2 e /22T scatter off photons in 
the thermal background and become e + e~ pairs, since the typical event rate for this process 
is much larger than that of other processes. On the other hand, for photons with energy 
smaller than the threshold energy ethr ~ m 2 /22T, this process is extremely suppressed by 
kinematics . As a result, radiative decay of X may cause photodissociation of light 

elements only if the decay of X occurs at T < m 2 e /22Q ) with Q being the threshold energy 
of photodissociation. For example, photodissociation of 4 He, whose binding energy is about 
20MeV, can be effective only for T < 500eV, i.e. for t > O(10 6 " 7 sec). 

The effects of the radiative decay of X essentially depend on the following parameters: 
the lifetime tx of X, the number density of X, and the energy e 7 o of the photons emitted in X 
decay. In this letter, for simplicity, we assume that X —>■ 7 + 7 with a 100% branching ratio, 
so e 7 o = mx/2. (For the effect of the hadronic decay of X, see Ref. [JO]].) To parameterize 
the abundance of X, we use the yield variable Yx, which is the ratio of the number density 
of X to that of photons: Yx = nx/n^. The yield variable is essentially the number of X in 
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a comoving volume, and it evolves with time as Yx = Yx§e~ t ' TX ■ Therefore, once Tx, Yxo, 
e 7 o, and rj are fixed, we can calculate the primordial abundances of the light elements. 

The procedure used in this letter is as follows. We first solve the Boltzmann equations 
for the distribution functions of the high energy photons and electrons in order to determine 
the abundance of cascade photons. Then, combining the derived cascade spectrum with the 
modified Kawano code |12j in which effects of the photodissociation are taken into account, 



we calculate the primordial abundances of the light elements, y 2p , ys p , and Y p , where y 2p and 
?/3 P are the primordial number fraction of D and 3 He relative to hydrogen H, respectively, 
and Y p is the primordial mass fraction of 4 He. (For details, see Refs. [§,0.) 

These theoretical predictions should be compared with the constraints obtained from 
observations. First, let us discuss the constraints on D and 3 He. Since D is only destroyed 



after BBN, the mass fraction of D decreases with time, and hence 13 



2/2, ism < V2 P R X , (l) 
2/2© < V2pRx, (2) 

where Rx = ^H, P /^H,now, 2/2 is the number fraction of D relative to H, X H is the mass 
fraction of H, and the subscripts "ism" , "0" , and "now" denote the abundances in interstellar 
matter, in the solar system, and in the present universe, respectively. Furthermore, taking 
into account the chemical evolution of D and 3 He, we obtain a third constraint ||13|| : 

I ~V2p + f — - lj V3p\ 2/20 - ~^ 2 P^30 



+ (yip + 2/2 P 2/3pj Rx < o, (3) 

where ys is the number fraction of 3 He relative to H, and g% is the survival fraction of 3 He 



in stellar processes. Usually, g 3 is estimated to be 0.25 - 0.5 JTJ 



From observation, the present abundances of the light elements (with 1-a error) are 



known to be [[L5|] 



y 2 ,i8m = (1-6 ± 0.2) x 10" 5 , (4) 



y 2Q = (2.57 ±0.92) x 10~ 5 , (5) 
2/3o = (1-52 ± 0.34) xl(T 5 , (6) 
R x = 1.1 ±0.04. (7) 

Here, we briefly comment on the detection of D in the primordial Hi cloud. The absorbed 
line observed in Ref. |Tj| suggests y% — 2.5 x 10 -4 , which is much larger than the value 



given in (|j) and (|||). However, it is claimed that the observed absorption line may be due 



to a Doppler-shifted hydrogen line |17j. Therefore, in this paper, we do not adopt the D 



abundance obtained in Ref. [16 



For fixed values of y 2p , y 3p , and g 3 , by regarding j/ 2 ,ism, 2/20, 2/30, and R x as statistical 
variables which obey Gaussian distributions, we can calculate the probability P(y2 P , 2/3 P ; #3) 
that the conditions flU) - @ are satisfied simultaneously: 



P(y2 P , 2/3 P ; 93) = / dy 2 ,ismdy2Qdy3 Q dRx 
x/(y 2 ,i Sm ; 2/2,ism, ^l 2M Jf(y2&; 2/20, crj 20 ) 

xfiVse; 27 3 0, a 2 y: Jf{R x ; R x , a 2 Rx ). (8) 

In Eq.(|8J), the integration is performed in the volume V in which all the constraints (|I|) - 
(§) are satisfied, and f(x; x, a 2 ) is the Gaussian distribution function for x with mean x and 
variance a 2 . 

By using the probability P(j/2p, yap] 9z), we can obtain a constraint on the primordial 
abundance of D and 3 He. For a fixed value of g 3 , we exclude the parameter set (2/2 P , 2/3p) if 
P(y2p, y3p] 93) < 0.05. Here, we note that our P(j/ 2 p, J/3 P ; #3) = 0.05 contour with g 3 = 0.25 
is almost the same as the 95% C.L. constraint of Hata et al. in Ref. |R|. That is, comparing 
the constraints on y 3p for a fixed value of j/ 2p , we checked that the discrepancy between the 
two approaches is at most (10 — 20)%; hence, we conclude that they are consistent. 

Next, we quote the constraint on the primordial abundance of 4 He, as estimated from 



the low-metallicity Hn regions 



Y p = 0.232 ± 0.003(stat) ± 0.005(syst). (9) 



We are now ready to show our numerical results. In our analysis, we compare the 
constraint on the baryon-to-photon ratio 77 derived from ([!]) - (|3]), with the contours which 
yield the correct primordial abundance of 4 He. In Fig. [I], we show the contour for the 5% 
probability constraints in the mxYxo vs. rj plane for several values of tx- Here, we take the 
mass of X to be ITeV, a neutron lifetime of 887.0sec |01, and g 3 = 0.25. Also, we show the 



contours which will yield specified abundances of 4 He. 

The typical behavior of the constraint can be understood as follows. If e^Yx is too small, 
then the abundance of cascade photons is so suppressed that the standard BBN scenario is 
almost unchanged. In this case, constraints from D and 3 He prefer 77 to be (3 — 7) x 10~ 10 , 
which is too large to be consistent with the 4 He constraint. However, if e^Yx has the correct 
value, then sufficient amounts of D and 3 He are destroyed by cascade photons, and a smaller 
value of 77 is allowed. (Remember that, in the standard BBN scenario, smaller values of 77 
give larger values of ?/2 P and ys p .) In particular, as can be seen in Fig. [I], Y p = 0.232 may be 
consistent with the constraints from D and 3 He for g 3 = 0.25. Furthermore, we also checked 
that, even for g 3 = 0.5, Y p = 0.235 is acceptable if e 7 o ~ lOMeV. Thus, if there is an exotic 
particle, we can have a better agreement between observation and the theoretical prediction. 
If e^Yxo is too large, then too much D and 3 He are destroyed; hence such a parameter region 
is excluded. Notice that the distribution function of cascade photons depends primarily on 
the total amount of injected energy || if the mass of X is much larger than e t h r - In this 
case, our results are almost independent of m x for fixed mxY X o- 

One may worry that too much 7 Li may be photodissociated in this scenario, since the 
threshold for the photodissociation of 7 Li is very low (~ 2.5MeV). However, the cross section 
for the photodissociation of 7 Li is about four times smaller than that of D dissociation, so 
the abundance of 7 Li is not reduced much compared with the D abundance. Numerically, 
we checked that (30-40)% of 7 Li is photodissociated in our scenario. For example, for 
?7 = 2 x 10~ 10 , we get n7 Li /nn = (0.7 — 1.0) x 10~ 10 , and a larger abundance of 7 Li is obtained 
if we use a larger or smaller value of 77. (Notice that we prefer 77 < 3 x 10~ 10 . For this range, 
7 Li production though 7 Be is not important. For rj > 3 x 10 -10 , we cannot predict 7 Li 



abundance since the cross section for 7 Be photofission is not available.) Taking into account 
the uncertainty in the observed abundance of 7 Li, viz. n7 Li /nn = (1.4^7) x 10~ 10 [20"|, the 



7 Li abundance in our scenario is consistent with observation. Furthermore, in our scenario 
6 Li can be produced via 7 Li(7, n) 6 Li. This reaction may produce much more 6 Li than 
standard BBN does. Numerically, we expect a few % of 7 Li can be converted into 6 Li for 
r x ^ 10 6 sec. This is consistent with present observation |2l| (n6 Li /n7 Li < 0.1), and may be 
checked in future observations. 

Let us now discuss the constraints on rx- First of all, the blackbody spectrum of CMBR 



observed by COBE [22] gives us a severe constraint on particles with lifetime longer than 



10 6 sec p3| , which is when the double Compton process (7 + e ^7 + 7 + e) freezes 



out. In particular, if e 70 Yxo ~ 10 9 GeV, then r x > 10 6 sec is forbidden [p3| . Furthermore, 



if Tx ^ 10 6 sec and e 70 ^ 20GeV, then cascade photons destroy 4 He effectively, which 
may result in the overproduction of D and 3 He. Thus, we have an upper bound on tx of 
~ 10 6 sec. On the other hand, as Tx decreases, the threshold for double photon pair creation 
e t hr becomes low, and the number of photons contributing to D and 3 He photodissociation 
decreases, for a fixed value of e l0 Y X o. In this case, a larger value of mxYxo is required in 
order to destroy sufficient amounts of D and 3 He, as can be seen in Fig. [I]. However, a 
larger initial mass density of X would speed up the expansion of the universe, so 4 He may 
be overproduced. These arguments exclude Tx shorter than ~ 10 4 sec. As a result, to make 
our scenario viable, we must adopt 10 4 sec < r x < 10 6 sec. In Fig. 0, we show the allowed 
region (i.e., the region where P(y2p, V3p', #3) > 0.05) on the t x vs. m x Yxo plane. Here, we 
use g 3 = 0.25, and we vary 77 so that Y p = 0.232 (dotted contour) or 0.235 (solid contour). 

Finally, we suggest several candidates for the unstable particle X, specifically within 
supergravity models. In such models, there may be several particles which have a long 
lifetime (10 4 sec ^ rx ~ 10 6 sec). Probably the most famous such particle is the gravitino, 
whose cosmological implications have been well-investigated . As we will see below, the 
gravitino may have the properties required to solve the crisis in BBN, although in the past it 
has been regarded as a source of cosmological difficulties (the so-called gravitino problem). 
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We can quantitatively check whether the gravitino can solve the difficulty in BBN, since its 
interactions are almost unambiguously determined p^| . If the gravitino decays only into a 



photon and photino pair, then its lifetime is estimated as 

t x ~ 4 x 10 5 sec x (lTeV/m 3/2 ) 3 , (10) 

where m 3 / 2 is the gravitino mass. Thus, if the gravitino mass is about ITeV, its lifetime 
is about 10 5_6 sec, which is appropriate for our purpose. Notice that the photino produced 
by the gravitino decay interacts very weakly with the thermal background. Thus, our 
results can be applied to the gravitino if we rescale the m^V^o-axis by a factor of ~ 0.5. 
Furthermore, assuming an inflationary universe, the gravitino abundance is determined just 
after the reheating period, and depends only on the reheat temperature T R ||: 

Fxo^2xlO- 11 (T i? /10 10 GeV). (11) 

For m 3/2 ~ ITeV, T R ~ 10 8 ~ 9 GeV gives e j0 Y X0 ~ 10-( 9 - 10 )GeV. We note here that 
Tr ~ 10 9 GeV may be realized in the chaotic inflation model [S5| if the inflaton field decays 
though gravitational interactions. 

If the gravitino is the lightest superparticle, then we can construct another scenario. In 
this case, the next-to-the lightest superparticle (NLSP), which we assume to be the photino, 
decays into a gravitino and a photon. By a loose tuning of the parameters, the NLSP 
can be identified with an unstable particle X which solves the difficulty in BBN. Thus, in 
supergravity models, we have several candidates for the unstable particle X which can make 
the BBN scenario viable. 

In summary, an exotic particle with lifetime tx ~ 10 4 ~ 6 sec can solve the crisis in BBN 
if mxYxo is tuned within a factor of 2-3. Candidates for X include the gravitino and the 
photino, both of which naturally appear in supergravity models. Of course, even in other 
types of models, one may be able to find a candidate for the radiatively decaying particle 
X. 
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FIG. 1. Constraints from 3 He and D, P(y2 P , ys P ', 53) > 0.05, shown in the 77 vs. mxYxo plane 
for mx = ITeV, and 53 = 0.25. The lifetime of X is take to be 3 x 10 3 sec (dotted line), 10 4 sec 
(long-dashed line), 10 5 sec (short-dashed line), and 10 6 sec (dotted-dashed line). The allowed regions 
are the interiors of the curves. The solid lines are the contours which yield Y p = 0.232, 0.235, and 
0.240 (from left to right) for t x = 3 x 10 3 sec. 
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FIG. 2. The contour for P(y2 P , U3 P ', 93) = 0.05 with mx = ITeV and gs = 0.25. rj is determined 
so that Yp = 0.232 (dotted) and 0.235 (solid). Notice that the interiors of the contours are favored. 
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